Two-photon probe of the Jaynes-Cummings 
model and symmetry breaking in circuit QED 
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Superconducting qubits 1 ' 2 behave as artificial two-level atoms and are used to investi- 
gate fundamental quantum phenomena. In this context, the study of multi-photon excita- 
tions-2^«ii£»£ occupies a central role. Moreover, coupling superconducting qubits to on-chip 
microwave resonators has given rise to the field of circuit QED^ ' 10 ' n i 12 i 13 i 14 i 15 . In con- 
trast to quantum-optical cavity QED 16 i 17 i 18 ] 19 ; circuit QED offers the tunability inherent 
to solid-state circuits. In this work, we report on the observation of key signatures of 
a two-photon driven Jaynes-Cummings model, which unveils the upconversion dynamics 
of a superconducting flux qubit 2 " coupled to an on-chip resonator. Our experiment and 
theoretical analysis show clear evidence for the coexistence of one- and two-photon driven 
level anticrossings of the qubit-resonator system. This results from the symmetry break- 
ing of the system Hamiltonian, when parity becomes a not well-defined property 21 . Our 
study provides deep insight into the interplay of multiphoton processes and symmetries in 
a qubit-resonator system. 

In cavity QED, a two-level atom interacts with the quantized modes of an optical or microwave cavity. 
The information on the coupled system is encoded both in the atom and in the cavity states. The latter 
can be accessed spectroscopically by measuring the transmission properties of the cavity—, whereas the 
former can be read out by suitable detector s 18 ! 19 . In circuit QED, the solid-state counterpart of cavity 
QED, the first category of experiments was implemented by measuring the microwave radiation emitted 
by a resonator (acting as a cavity) strongly coupled to a charge qubit^. In a dual experiment, the state of 
a flux qubit was detected with a DC superconducting quantum interference device (SQUID) and vacuum 
Rabi oscillations were observed—. More recently, both approaches have been exploited to extend the 
toolbox of quantum optics on a chip i 1 ? 12 ! 13 ! 14 ! 15 ! 22 . Whereas all these experiments employ one-photon 
driving of the coupled qubit-resonator system, multi-photon studies are available only for sideband tran- 
sitions-^ or bare qubits^^^. The experiments discussed in this work explore, to our knowledge for the 
first time, the physics of the two-photon driven Jaynes-Cummings dynamics in circuit QED. In this con- 
text, we show that the dispersive interaction between the qubit and the two-photon driving enables real 
level transitions. The nature of our experiment can be understood as an upconversion mechanism, which 
transforms the two-photon coherent driving into single photons of the Jaynes-Cummings dynamics. This 
process requires energy conservation and a not well-defined parity— of the interaction Hamiltonian due 
to the symmetry breaking of the qubit potential. Our experimental findings reveal that such symme- 
try breaking can be obtained either by choosing a suitable qubit operation point or by the presence of 
additional spurious fluctuators^. 

The main elements of our setup, shown in Figs.[l^ and b, are a three- Josephson-junction flux qubit, 
an LC-resonator, a DC SQUID and a microwave antenn a 24 ! 25 . The qubit is operated near the optimal 
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Figure 1: Experimental architecture and theoretical model, a, Schematic representation of the circuit. The 
Josephson junctions are represented by crosses. The readout DC SQUID (black rectangle with two crosses), which is 
coupled to the three- Joscphson-j unction flux qubit (red) by a purely geometric mutual inductance, is shunted with an 
LC circuit acting as a quantized resonator (blue). All circuit elements within the shaded area are at the base temperature 
T ~ 50 mK of a dilution refrigerator. Microwave signals can be applied to qubit and resonator via an on-chip antenna 
(green). The signal-to- noise ratio of the driving is increased by means of cold attenuators, b, Scanning electron microscopy 
micrograph of the flux qubit and the readout DC SQUID, c, Pulse scheme for qubit microwave spectroscopy with the 
adiabatic shift pulse method. First, the qubit is biased at a suitable readout point (<J> X ^ 1.5$o) using a superconducting 
coil. Then, it is initialized in its ground state |g) by waiting a sufficiently long time. In the next step, we use a rectangular 
pulse ( "adiab. shift pulse" , green) to shift the qubit adiabatically to the desired operation point. There, it is irradiated with 
a 100ns spectroscopy pulse ("microwave pulse", green) of frequency u)/2-k. Both shift and spectroscopy pulse are applied 
via the microwave antenna. After the end of the shift pulse, a readout pulse to the DC SQUID measurement lines probes 
the qubit state. By averaging over thousands of such measurements, the probability P a to find the qubit in the excited state 
c) is calculated. Using the readout protocol described above, the qubit state can be detected via the DC SQUID switching 
signal also at the optimal point, despite the fact that there the mean value I p (& z ) of the circulating current in the qubit 
loop vanishes^, d, Sketch of the upconversion dynamics describing the physics that governs the experiments discussed 
in this work, see Eq. JT]|. The qubit (red) has a level splitting fcw q and the resonator (blue) frequency is w T /2ir. We only 
consider the relevant case of two-photon driving (green) and assume that the system predominantly decays via the resonator. 
The two-photon driving frequency is w and the qubit-resonator coupling strength is gsinS, where sinf? = A/o; r ~ 0.63. 
Depending on <t> x , the qubit potential (red double well; the x-axis represents the phase variable ipm) is either symmetric 
or does not have a well-defined symmetry. Consequently, level transitions are allowed or forbidden respectively. 



flux bias <& x = 1.5i»o and can be described with the Hamiltonian H q = (ea z + Aa x ) /2, where a x and a z 
are Pauli operators. From low-level microwave spectroscopy we estimate a qubit gap A/h = 3.89 GHz. 
By changing $ x , the quantity e = 2I p (<I> X — 1.5$n) and, in turn, the level splitting k\o q = Ve 2 + A 2 can 
be controlled. Here, ±/ p are the clockwise and counterclockwise circulating persistent currents associated 
with the eigenstates |±) of ea z . Far away from the optimal point, |±) correspond to the eigenstates |g) and 
|e) of Hq. The qubit is inductively coupled to a lumped-element LC-resonator, which can be represented 
by a quantum harmonic oscillator, H T — hw r {a^a + 1/2), with photon number states |0), |1), |2) , . . . and 
boson creation and annihilation operators d^ and d respectively This resonator is designed such that its 
fundamental frequency, cj r /27r = 6.16 GHz, is largely detuned from A/h. The qubit-resonator interaction 
Hamiltonian is H q Y = hga z (d^ + d), where g = 2ir x 115 MHz is the coupling strength. The iC-circuit 
also constitutes a crucial part of the electromagnetic environment of the readout DC SQUID. In this 
way, the flux signal associated with the qubit states |±) can be detected while maintaining reasonable 
coherence times and measurement fidelit y 24 i 25 . 

In order to probe the properties of our system, we perform qubit microwave spectroscopy using an 
adiabatic shift pulse techniqu e 24 ' 25 (Fig.[TJ;). The main results are shown in Figs. [5^, and b. First, there 
is a flux-independent feature at approximately 6 GHz due to the resonator. Second, we observe two 
hyperbolas with minima near 4 GHz ~ A/h and 2 GHz ~ A/2h, one with a broad and the other with 
a narrow lincwidth. They correspond to the one-photon (uj = uj q ) and two-photon (2uj = uj q ) resonance 
condition between the qubit and the external microwave field. Additionally, the signatures of two-photon 
driven blue sideband transitions are partially visible. One can be attributed to the resonator, \g,0) — > 
|e, 1), and the other to a spurious fluctuator— . We assume that the latter is represented by the flux- 
independent Hamiltonian H[ = (e*&* + A* a*) /2 and coupled to the qubit via H q i = hg*a z a z , where 
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Figure 2: Qubit microwave spectroscopy: Data and simulations, a, Center frequency of the measured absorption 
peaks (symbols) plotted versus the flux bias: Red dots - qubit one-photon signal; green triangles - qubit two-photon signal; 
grey open triangles — resonator; blue squares - qubit-resonator blue sideband signal; blue diamonds - qubit-fluctuator blue 
sideband signal. The lines are fits to the data based on the undriven Hamiltonian H u . Of particular interest are the 
large and the small anticrossing for to ss lu t and 2lu fa uj r respectively. This constitutes direct evidence that the two-photon 
driving selectively excites only the qubit, but is strongly suppressed for the cavity. Consequently, the vacuum Rabi coupling 
g can be obtained from two-photon spectroscopy. On the contrary, the one-photon driving populates the cavity and gives 

rise to the enhanced coupling gy (N). b, The measured probability P c to find the qubit in the excited state plotted as 
a function of the flux bias and the microwave excitation frequency. The black box denotes the area shown in Fig. [3^,. c, 
Probability P e for the driven system described by obtained from numerical simulations using the time-trace-averaging 
method. The parameters are derived from the fit in a. In our simulations, starting from the ground state |g), the average 
over a full 100 ns time trace consisting of 10000 discrete time points already gives excellent agreement with the experimental 
data of b (see main text). Since there are no terms describing dissipation in H^, the linewidth of the peaks is caused by 
power broadening. The values of the mutual inductances are based on numerical estimates. The black box denotes the 
area shown in Fig. 0] 



<7* and a* are Pauli operators. Exploiting the different response of the system in the anticrossing region 
under one- and two-photon driving, as explained in Fig. [2^, the center frequencies of the spectroscopic 
peaks can be accurately fitted to the undriven Hamiltonian H u = H q +H T +H{+H q r +H q ,{. Setting e* = 

(cf. methods), we obtain g/2n = 115MHz, (N) ~ 10, I p = 367 nA, w f /27r = ^e* 2 + A* 2 /h = 3.94 GHz 
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Figure 3: Qubit microwave spectroscopy close to the qubit-resonator anticrossing under two-photon driv- 
ing, a, Measured probability P c to find the qubit in the excited state plotted vs. the applied flux and the microwave 
excitation frequency. The black rectangle denotes the area shown in c and e. The solid lines represent the fit obtained by 
the undriven Hamiltonian H u . b, Maximum height of the spectroscopy peaks under one- photon (red dots) and two-photon 
(green triangles) driving plotted vs. the flux bias. The solid lines are guides to the eye. c, Probability P e obtained from 
simulations including driving, but no dissipation (time-trace-averaging method). The parameters are derived from the nu- 
merical fit of Fig. The signature of an anticrossing is visible, d, The green curve shows the split-peak profile of P c along 
the vertical line in c. The blue line shows the single-peak result obtained for the same flux bias, when turning off the qubit- 
resonator coupling by setting g = 0. e, Probability P c obtained from simulations including driving and dissipation based 
on the Lindblad formalism^. For simplicity, the spurious fluctuator is not included here. We assume a qubit relaxation 
time Ti = 300 ns, a dephasing time T v = f5ns and a resonator quality factor Q = lu t /k = 2tt X 6. f6 GHz/400 MHz ~ fOO 
corresponding to a resonator decay time r ~ 3 ns. The spectroscopy signal is chosen to be the average over the last 20 ns 
of a 100 ns time trace. When the qubit and the resonator become degenerate, the spectroscopy signal fades away, f, Again, 
the green curve shows the split-peak profile of P c along the vertical line in e. The blue line shows the single-peak result 
obtained for the same flux bias when setting g = 0. Differently from the nondissipative case (c and d), the single peak is 
approximately ten times higher than either one of the strongly reduced split peaks (green curve). This demonstrates that 
the vanishing of the spectroscopy signal observed in the experimental data (cf. a, b) and in the simulations with dissipation 
(cf. e) is not caused by qubit decoherence. 



and 5*sin6»* = 37 MHz, where sin6>* = A*/hu f . 

Further insight into our experimental results can be gained by numerical spectroscopy simulations 
based on the driven Hamiltonian H,i = H u + ff m ,q + #m,r + H 1TL j. Here, H n i,q = ^<J z cosujt, H m ,x = 
rj (a) + a) cos cot and H m ,{ = ^-a* costut represent the driving of the qubit, resonator and fluctuator 
respectively. We approximate the steady state with the time average of the probability P c to find the 
qubit in |c) (time-trace-averaging method). Inspecting Fig. we find that for the driving strengths 
il/h = 244 MHz, rj/h = 655 MHz and O* = our simulations match well all the experimental features 
discussed above. Using rj and the relation (TV) = {t)/k) for the steady-state mean number of photons 
of a driven dissipative cavity, we estimate a cavity decay rate k ~ 210 MHz. This result is of the same 
order as k ~ 400 MHz estimated directly from the experimental linewidth of the resonator peak. The 
large k is due to the galvanic connection of the resonator to the DC SQUID measurement lines (Fig.[TJi). 

To elucidate the two-photon driving physics of the qubit-resonator system we consider the spec- 
troscopy data near the corresponding anticrossing shown in Fig. For 2lo = uj q = u> r , the split 
peaks cannot be observed directly because the spectroscopy signal is decreased below the noise floor 
SP C ~ 1 — 2%. This results from the fact that the resonator cannot absorb a two-photon driving and 
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its excitation energy is rapidly lost to the environment (k > g/2n). In contrast, for the one-photon 
case (oj = Lo q — u> T ), there is a driving-induced steady-state population of (N) ~ 10 photons in the 
cavity. Accordingly, the one-photon peak height shows a reduction by a factor of approximately two, 
whereas the two-photon peak almost vanishes, see Fig. [3h>. To prove that this effect is only due to the 
resonator, we compare the simulation results from the time-trace-averaging method to those obtained 
with the standard Lindblad dissipative bath approach (Figs. [5J>f). Altogether, our experimental data 
and numerical simulations constitute clear evidence for the presence of a qubit-rcsonator anticrossing 
under two-photon driving. 

The second-order effective Hamiltonian under two-photon driving can be derived using a Dyson-series 
approach (cf. methods). Starting from the first-order driven Hamiltonian Ha and neglecting the cavity 
driving and the fluctuator because of large-detuning conditions, we obtain 

H m = ^^ + ^ 8 m2flcosfl(a + e- <awt + fr_e +awt ) 

— Kg sin o(& + a + &_a)^j + Hu> x (j^® > W 

where a + and a_ are the qubit raising and lowering operators, sin 9 = A/u q and cos 9 = e/w q . The 
upconversion dynamics sketched in Fig.[TJi is clearly described by Eq. |T|). The first two terms represent 
the qubit and its coherent two-photon driving with angular frequency u. The last two terms show 
the population transfer via the Jaynes-Cummings interaction to the resonator. The Jaynes-Cummings 
interaction in this form is valid only near the anticrossings (9 ~ 9, 9 = ~ 0.63; cf. methods). As 
discussed before, the resonator will then decay emitting radiation of angular frequency 2oj. 

The model outlined above allows us to unveil the symmetry properties of our system. Even though 
the two-photon coherent driving is largely detuned, ^ = u 3> § sin#, a not well-defined symmetry of the 
qubit potential permits level transitions away from the optimal point. Because of energy conservation, 
i.e. frequency matching, these transitions are real and can be used to probe the qubit-resonator anti- 
crossing. The effective two- photon qubit driving strength, (O 2 sin 2 9/4 A) cos 9, has the typical structure 
of a second-order dispersive interaction with the extra factor cos 9. The latter causes this coupling to 
disappear at the optimal point. There, the qubit potential is symmetric and the parity of the interaction 
operator is well defined. Consequently, selection rules similar to those governing electric dipole transi- 
tions hold—. This is best understood in our analytical two-level model, where the first-order Hamiltonian 
for the driven diagonalized qubit becomes H Q p = T^z + T®x (e +luJt + e~ Mt ) at the optimal point. In 
this case, one-photon transitions are allowed because the driving couples to the qubit via the odd-parity 
operator a x . In contrast, the two-photon driving effectively couples via the second-order Hamiltonian 
Hop = T^z + T&K^z (c +lu)t + e~ lwt ) . Since a z is an even-parity operator, real level transitions are 
forbidden (cf. methods). We note that the second <r z -term of H Q p renormalizes the qubit transition fre- 
quency slightly and can be neglected in Eq. ([T]), which describes the real level transitions corresponding 
to our spectroscopy peaks. The intimate nature of the symmetry breaking resides in the coexistence 
of a x - and o^-operators in the first-order Hamiltonian H^, which produces a nonvanishing o^-term in 
the second-order Hamiltonian of Eq. (p}. This scenario can also be realized at the qubit optimal 
point by the fluctuator terms a* and a*. As illustrated in Fig. [U their presence causes a revival of the 
two-photon signal and the discussed strict selection rules no longer apply. Accordingly, we observe only 
a reduction instead of a complete suppression of the two-photon peaks near the qubit optimal point in 
the experimental data of Fig. \Bp. 

In conclusion, we use two-photon qubit spectroscopy to study the interaction of a superconducting 
flux qubit with an LC-resonator. We show experimental evidence for the presence of an anticrossing 
under two-photon driving, permitting us to estimate the vacuum Rabi coupling. Our experiments and 
theoretical analysis shed new light on the fundamental symmetry properties of quantum circuits and the 
nonlinear dynamics inherent to circuit QED. This can be exploited in a wide range of applications such 
as parametric up-conversion, generation of microwave single photons on deman d 11 i 26 i 27 or squeezing—. 

APPENDIX 



Two-photon driven Jaynes-Cummings model via Dyson series. We now derive the effective second-order 
Hamiltonian describing the physics relevant for the analysis of the two-photon driven system. We start 
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Figure 4: Two-photon spectroscopy simulations close to the optimal point using the time-trace-averaging 
method, a, Probability P e plotted vs. driving frequency and flux bias. The parameters are the same as those used 
in Fig. In particular, the fluctuator parameters are e* = <-> sin#* = 1 and Q* = 0. The spectroscopy signal 
vanishes completely at the optimal point (<E> X = LS^o) because of the specific selection rules associated with the symmetry 
properties of the Hamiltonian— L. b, Same as in a, however, for sin 8* = 0.3 and Q* = 280 MHz. Here, the coexistence of 
the flux-independent first-order &*- and cr*-terms of the fluctuator gives rise to a nonvanishing second-order crj-term, even 
at the qubit optimal point. In other words, the presence of the fluctuator breaks the symmetry of the total system at the 
optimal point and, again, parity becomes not well defined. Consequently, the spectroscopy signal is partially revived. In 
reality, there is a high probability that not a single, but an ensemble of fluctuators with some distribution of frequencies 
and coupling strengths contributes to the symmetry breaking. Furthermore, when the experimental resolution is limited, 
a single peak will be detected and detailed structure such as the one of Fig|4p will not be visible. This is actually the case 
in our measurements (Fig. [2p). 



from the first-order Hamiltonian in the basis |±) 



e A / + 1\ 

H = 2 (Xz + ~2 <Tx+ T \ + 2) 



+ hg a z (a^ + a) + —a z cos cut . (2) 

Here, in comparison to Hd, the terms associated with the fluctuator are not included (e* = A* = SI* = 0) 
because the important features are contained in the driven qubit-resonator system. Additionally, we focus 
on the two- photon resonance condition uj t — uj q = 2uj. Thus, the driving angular frequency ui is largely 
detuned from uj r and the corresponding term in can be neglected (rj = 0). Next, we transform the 
qubit into its energy eigenframe and move to the interaction picture with respect to qubit and resonator, 
<7± — > <T±e ±I " qt , a — ► ae _I " r * and a) —> a) 'e +luJlt . After a rotating wave approximation we identify 
the expression S^e +UJt + Sc~ luJt , where the superoperator S = § (cos 9a z — singer.) and its Hermitian 
conjugate = ^ (cos(9<r z — sin#er + ). In our experiments the two-photon driving of the qubit is weak, 
i.e. the large-detuning condition uj q — oj = oj ^> ^ sin 9 is fulfilled. In such a situation, it can be 
shown that the Dyson series for the evolution operator associated with the time-dependent Hamiltonian 
—hgsm8(a + a + a_a)) + (Sc^' lult + S^e +iut ) can be rewritten in an exponential form U = e -1 ^*''/^, 
where 



H cS = - hg sind (& + ae +lSt + a_a)e- lSt 

= - hg sin9 (a + ac +m + c^aV* 5 *) 

Q, 2 / 1 \ 

+ — f sin 2 6 cos 6 a x + - sin 3 9a z j. (3) 

Here, 6 = cj q — uj t is the qubit-resonator detuning. In Eq. |(3J|, the dispersive shift ^ sin 3 9a z is a reminis- 
cence of the full second-order <r z -component of the interaction Hamiltonian, sin 3 9 a z (e+ lut + e~ 4Wt ) 2 . 
The terms proportional to a z exp ±l2uJt are neglected implicitly by a rotating wave approximation when 
deriving the effective Hamiltonian H e g of Eq. In this equation, the <r z -term renormalizes the qubit 
transition frequency, and, in the vicinity of the anticrossing (\S\ < gsin9, sin9 = — ~ 0.63), the Hamil- 
tonian 

of Eq. (P) can be considered equivalent to H eS . In this situation, the symmetries of the 
system are broken and our experiments demonstrate the existence of real level transitions. 

Selection rules. The potential of the three- Josephson-junction flux qubit can be reduced to a one- 
dimensional double well with respect to the phase variable ip m — . At the optimal point ($ x = 1.5<I>o)j this 
potential is a symmetric function of (p m . For our experimental parameters, we can assume an effective 
two-level system. The two lowest energy eigenstatcs |g) and |e) are, respectively, symmetric and anti- 
symmetric superpositions of |+) and |— ). Thus, |g) has even parity and |e) is odd. In this situation, the 
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parity operator A can be denned via the relations II |g) = + |g) and ft |e) = — |c). The Hamiltonian of the 
classically driven qubit is —a z + § cosujta x . For a one-photon driving, uj = A/h (energy conservation), 
the Hamiltonian in the interaction picture is ^cr x , where a x = |g) (e| + |e) (g|. This is an odd-parity oper- 
ator because the anticommutator {IT, a x } = and, consequently, one-photon transitions are allowed. For 
a two- photon driving, uj = A/2h (energy conservation), the effective interaction Hamiltonian becomes 
!^<j 2 , where a z = |c) (c| — |g) (g|. Since the commutator [ft, & z ] = 0, this is an even-parity operator and 
two-photon transitions are forbidden* 2 ^. These selection rules are analogous to those governing electric 
dipolc transitions in quantum optics. On the contrary, in circuit QED the qubit can be biased away 
from some optimal point. In this case, the symmetry is broken and the discussed selection rules do not 
hold. Instead, we find the finite transition matrix elements ^ sin 9 and sin 2 9 cos 9 for the one- and 
two-photon process respectively. Beyond the two-level approximation, the selection rules for a flux qubit 
at the optimal point are best understood by the observation that the double-well potential is symmetric 
there (Fig. QJ1). Hence, the interaction operator of the one-photon driving is odd with respect to the 
phase variable tp m of the qubit potentia l 20 i 21 , whereas the one of the two-photon driving is even. Away 
from the optimal point ($ x ^ l.b^o), the qubit potential has no well-defined symmetry and no selection 
rules apply. 

Spurious fluctuators. The presence of spurious fluctuators in qubits based on Josephson junctions 
has already been reported previously—. In principle, such fluctuators can be either resonators or two- 
level systems. Since our experimental data does not allow us to distinguish between these two cases, for 
simplicity, we assume a two-level system in the simulations. In the numerical fit shown in Fig. [5^, we 
choose e* = due to the limited experimental resolution. Consequently, the coupling constant estimated 
from the undriven fit is not g*, but g*sm9*. Away from the qubit optimal point, especially near the 
qubit-resonator anticrossings, the effect of the observed fluctuator can be neglected within the scope of 
this study. Near the optimal point, its effect on the symmetry properties of the system can be explained 
following similar arguments as given above for the flux qubit. However, it is important to note that, 
differently from sin^ and cos 9, the fluctuator parameters sin#* and cos 9* are constants, i.e. they do not 
depend on the quasi-static flux bias $ x . 
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